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HIGH-RESOLUTION ASYMPTOTICS FOR THE ANGULAR 
BISPECTRUM OF SPHERICAL RANDOM FIELDS 1 

By Domenico Marinucci 

Universita di Roma "Tor Vergata" 

In this paper we study the asymptotic behavior of the angular 
bispectrum of spherical random fields. Here, the asymptotic theory is 
developed in the framework of fixed-radius fields, which are observed 
with increasing resolution as the sample size grows. The results we 
present are then exploited in a set of procedures aimed at testing non- 
Gaussianity; for these statistics, we are able to show convergence to 
functionals of standard Brownian motion under the null hypothesis. 
Analytic results are also presented on the behavior of the tests in the 
presence of a broad class of non-Gaussian alternatives. The issue of 
testing for non-Gaussianity on spherical random fields has recently 
gained enormous empirical importance, especially in connection with 
the statistical analysis of cosmic microwave background radiation. 

1. Introduction. Several statistical challenges are now arising in connec- 
tion with cosmological data, and more precisely, for the analysis of cosmic 
microwave background radiation (CMB). CMB can be viewed as a snapshot 
of the Universe approximately 3 x 10 5 years after the Big Bang; techno- 
logical progress has made possible a number of experiments aimed at mea- 
suring the properties of this radiation. Pioneering results were released in 
1992 by the NASA mission COBE [31], which was the first to release a full- 
sky map of CMB fluctuations; the statistical properties of these fluctuations 
were then further investigated by several balloon experiments, starting with 
BOOMERanG [8] and MAXIMA [14]. A major breakthrough is associated 
with two satellite missions, namely WMAP [4], which released its first data 
set in February 2003, with much more detailed data to come in the next 
four years, and Planck, which is scheduled to be launched in Spring 2007 
and expected to provide maps with much greater resolution. Over the next 
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ten years, an immense amount of cosmological information is expected from 
these huge data sets; statistical efforts needed to extract this information 
are equally challenging and impressive. 

From the mathematical point of view, CMB can be represented as a ran- 
dom field T(6,ip) indexed by the unit sphere S 2 , that is, for each azimuth 
< < 7r and elongation < (p < 2ir, T(9,ip) is a real random variable de- 
fined on some probability space. We shall always assume that T(9, (p) is 
a zero-mean, finite-variance, mean-square continuous and isotropic random 
field, that is, its distribution is invariant with respect to the group of rota- 
tions. Until very recently, the assumption of an isotropic Universe has been 
taken for granted in cosmological physics, as a consequence of Einstein's 
Cosmological Principle that the Universe should "appear the same" to an 
observer located anywhere in space. Quite interestingly, the first release of 
WMAP has raised some doubts on this condition [10, 15, 27]. Testing for 
isotropy is a very interesting topic, almost completely open for statistical 
research; we consider this issue, however, beyond the purpose of the present 
work. 

If an isotropic field is Gaussian, its dependence structure is completely 
identified by the angular correlation function and its harmonic transform, 
the angular power spectrum (to be defined in the next section). For non- 
Gaussian fields, the dependence structure becomes much richer, and higher- 
order correlation functions are of interest; in turn, this leads to the analysis of 
so-called higher-order angular power spectra. Because these angular power 
spectra are identically zero for Gaussian fields, they also provide natural 
tools to test for non-Gaussianity: this is a topic of greatest importance in 
modern cosmological data analysis. Indeed, on one hand the validation of 
the Gaussian assumption is urged by the necessity to provide a firm basis 
for statistical inference on cosmological parameters, which is dominated by 
likelihood approaches. More importantly, tests for Gaussianity are needed 
to discriminate among competing scenarios for the physics of the primordial 
epochs: here, the currently favored inflationary models predict (very close to) 
Gaussian CMB fluctuations, whereas other models yield different observa- 
tional consequences [3, 28, 29]. Tests for non-Gaussianity are also powerful 
tools to detect systematic effects in the outcome of the experiments. For 
these reasons, many papers have focused on testing for non-Gaussianity on 
CMB, some of them by means of topological properties of Gaussian fields 
(e.g., [9, 13, 26, 36, 37]), others through spherical wavelets [2], and still oth- 
ers by harmonic space methods (e.g., [16, 20, 21, 24]). A short survey of the 
literature on testing for non-Gaussianity on CMB is in [23]. 

In this paper, we investigate the asymptotic properties for the observed 
bispectrum of spherical Gaussian fields, and we analyze its use as a probe 
of non-Gaussian features. The bispectrum (defined in Section 2) is proba- 
bly the single most popular statistic to search for non-Gaussianity in CMB 
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data; on one hand, in fact, working on harmonic space is extremely con- 
venient, and the bispectrum is the simplest harmonic space statistic which 
is sensitive to non-Gaussian features. On the other hand, it is possible to 
derive analytically the behavior of the bispectrum for non-Gaussian fields of 
physical interest. In fact, although the procedures considered in the present 
work are new (to the best of our knowledge), a number of insightful and 
important papers have already considered the bispectrum for CMB data 
analysis, for instance, [19, 20, 21, 30]. Few analytical results, however, have 
so far been produced on the statistical properties of these procedures. The 
reason for this can be partially explained as follows: as described in the next 
section, the bispectrum is a function of some spherical harmonic coefficients 
o/ m 's; in the presence of an ideal experiment, the latter are easily derived 
from a map of CMB fluctuations by a harmonic transform performed on 
the observed data. However, it is important to stress that in realistic sit- 
uations the aim's are observed with error, due to instrumental noise, gaps 
in the maps and many other sources (these problems have been described 
in [33]). Determining the properties of the procedures in realistic situations 
that take into account all the features of real life experiments is extremely 
difficult; most work in CMB is thus based on comparisons of estimates based 
on real data with the expected results of simulations under a particular null 
hypothesis. In this paper, we assume that the a; m 's are observed without 
error; this is a simplifying assumption, which we adopt because it seems 
important to narrow the gap between data analysis practice and its math- 
ematical foundations, at least in an idealized case. Future work, however, 
should be directed at relaxing this assumption. 

In Section 2 we define the bispectrum, taking particular care to discuss 
the conditions to ensure that it represents a rotationally invariant statistic. 
The asymptotic behavior of its higher-order moments is considered in Sec- 
tion 3: these results are derived under Gaussianity, but the general technique 
by which they are established (which adopts a formalism from graph theory) 
may have some independent interest under broader assumptions. Section 4 
considers the effect of an unknown angular power spectrum, whereas Sec- 
tion 5 discusses statistical applications, with asymptotic results in Gaussian 
and non-Gaussian circumstances. The results we present in this section sug- 
gest that consistent tests of Gaussianity can exist even for random fields 
defined on a bounded domain, which is to some extent unexpected. Section 
6 discusses the rationale behind the approach presented and draws some 
conclusions; some technical results are collected in the Appendix. 

2. The angular bispectrum. The Fourier transform on the sphere is de- 
fined by the spherical harmonics, which can be written explicitly as 




for m > 
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Y lm (9, ip) = (-l) ro ^*_ m (0, tp) for m < 0, 

where the asterisk denotes complex conjugation and Pi m (cos9) denotes the 
associated Legendre polynomial of degree I, m, that is, 

im i f il 

p lm ( X ) = (-ir a - x^—p^x), nix) = W[ -^{* 2 - 1)', 

m = 0, 1, . . . , I, I = 1, 2,3, 

A detailed discussion of the properties of the spherical harmonics can be 
found in ([34], Chapter 5), or in [35]. For isotropic fields, the following spec- 
tral representation holds in the mean-square sense (see also [1, 22, 38]): 

oo I 

(1) T{6 ) v) = Y. E a lm Y lm (6,ip), 

1=1 m=—l 

where the triangular array {ai m } represents a set of random coefficients, 
which can be obtained from T(8, (p) through the inversion formula 

ai m = f W rT(e,p)YC m (e,p) S m9d6dp, 

(2) J -TV Jo 

m = 0,+l, ... ,+l,l = 1,2, ... . 

These coefficients are complex- valued, zero-mean and uncorrelated; hence, if 
T(9, ip) is Gaussian, they have a complex Gaussian distribution, and they are 
independent over I and m > [although a\_ m = (— l) m O;* m ], with variance 
E\ 

a im\ — Ci, m — 0, il, . . . , ±Z. The index / is usually labeled a multipole 
and in principle it runs from 1 to infinity; each multipole corresponds ap- 
proximately to an angular resolution of 180//°. In any realistic experiment, 
however, there is an upper limit (which we denote by L) on the multipoles 
we may observe, depending upon the resolution of the experiment and the 
presence of noise; L is reckoned to be of the order of 600/800 for WMAP 
and 2000/2500 for Planck. Strictly speaking, in CMB cosmology the spectral 
representation (1) is really only defined for I > 2; the so-called dipole / = 1 
is in fact dominated by kinematic effects and thus it is removed from the 
data. 

The sequence {C{\ denotes the angular power spectrum: we shall always 
assume that C\ is strictly positive, for all values of I. This condition is very 
mild; to draw an analogy with the theory of stationary random fields defined 
on R d , it is equivalent to the (very common) assumption that their spectral 
density is strictly positive at all frequencies. As discussed in Section 5, the 
condition is met by virtually all models of cosmological interest. A natural 
estimator for C\ is 

1 1 2 

(3) 21 + 1 ' a ' m ' ' ^ = 1,2,..., 
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which is clearly unbiased; see also [12]. As mentioned in the Introduction, 
if the field is Gaussian, the angular power spectrum completely identifies 
its dependence structure. For non-Gaussian fields, the dependence structure 
becomes much richer, and higher-order moments of the az m 's are of interest; 
this leads to the analysis of so-called higher-order angular power spectra. 

Generally speaking, the angular bispectrum can be viewed as the har- 
monic transform of the three-point angular correlation function, much as 
the angular power spectrum is the Legendre transform of the (two-point) 
angular correlation function. More precisely, write fij = (6i, <pi), for i = 1, 2, 3; 
we have 

sr(o 1 )r(o 2 )r(^ 3 ) 

oo 

(4) = E E B Z\T m3Y hmi WY hma (n 2 )Y l3m3 (n 3 ), 

Ji,I 2 ,i3=imi,m2,m3 
where the bispectrum B™i™ 2m3 is given by 

(5) ^hhl 3 2m3 = E( a hm 1 0'l 2 m2 a l3mi)- 

Here, and in the sequel, the sums over rrii run from — Zj to Zj, unless otherwise 
indicated. Both (4) and (5) are clearly equal to zero for zero-mean Gaussian 
fields. Moreover, the assumption that the CMB random field is statistically 
isotropic entails that the right- and left-hand sides of (4) should be left 
unaltered by a rotation of the coordinate system. Therefore I?™^™ 2 " 13 must 
take values ensuring that the three-point correlation function on the left- 
hand side of (4) remains unchanged if the three directions J2i,fi2 an d ^3 
are rotated by the same angle. Careful choices of the orientations entail that 
the angular bispectrum of an isotropic field can be nonzero only if: 

(a) h,h and Z3 satisfy the triangle rule, U < lj + Ik for all choices of 
hj-i k = 1) 2, 3, 

(b) h + h + h= even and 

(c) mi + ni2 + m 3 = 0. 

More generally, Hu [17] shows that a necessary and sufficient condition for 
B™fel 3 i2m3 to represent the angular bispectrum of an isotropic random field 
is that there exist a real symmetric function of Zi , Z2 , h , which we denote 
b^hhi such that we have the identity 

(6) KtZ 2m3 = ^hiZ 2m3b h^ 

bi x i 2 i z is labeled the reduced bispectrum. In (6) we are using the Gaunt 
integral Q^ 2m \ defined by 

^T™ = ££ n Y hmi (9^)Y l2m2 (9^)Y hm3 (e^)smed^M 

/ (2Z 1 + l)(2Z 2 + l)(2Z 3 + l) \ 1/2 /'Z 1 h h\(h h h\ 
\ 4vr J \0 0y v mi m 2 m 3 J ' 
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where the so-called "Wigner 3j symbols" appearing on the second line are 
defined in the Appendix. It can be shown that the Gaunt integral is iden- 
tically equal to zero unless the conditions (a)-(c) are fulfilled. Often the 
dependence on mi , 777.2 , 7773 , which does not carry any physical information 
if the field is isotropic, is eliminated by focusing on the angular averaged 
bispectrum, defined by 

g . \ " \ ~> \ "* f h h h \ ^m 1 m 2 ni3 

hhh 2^ 2^ \ mi m2 m3 I hhh 



(7) m 1 =-l 1 m 2 =-l2m 3 =-l, 

•(2Zi + l)(2/ 2 + l)(2Z 3 + l)\ 1/2 fh l 2 h 



4tt J \0 



where we have used (48) (see the Appendix). In practice, of course, the 
bispectrum is not observable; its minimum mean-square error estimator is 
provided by Hu [17], 



mi=—h rri2=—l2 m3=—h 



The statistic B^i^ is called the (sample) angle averaged bispectrum; for 
any realization of the random field T, it is a real- valued scalar, which does 
not depend on the choice of the coordinate axes, and it is invariant with 
respect to permutation of its arguments 

1/2 

Now note that, under Gaussianity, the distribution of ai m /C, does not 
depend on any nuisance parameter. The bispectrum can hence be easily 
made model-independent; namely, we can focus on the normalized bispec- 
trum, which we define by 

(8) I hhl3 = (-l)(^+'3)/2 g 'ifa'a 



V C h C h C h 

The factor (— l)('i+'2+'3)/2 j g usua Hy no t included in the definition of the 
normalized bispectrum; it corresponds, however, to the sign of Wigner's 
coefficients for mi = 7712 = 777.3 = 0, and thus it seems natural to include it to 
ensure that ^ 1 j 2 j 3 and 6^/3/3 share the same parity [see (7)]. 

In practice, of course, ^ 1 / 2 / 3 is infeasible, and it is thus replaced by the 
statistic 



\ hh =(-i) {h+h+h)/2 ! BlM \ 

/C h Ci 2 Ci 3 



see [12]. 
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3. Higher-order moments of the angular bispectrum. In this section we 
shall investigate the behavior of the higher-order moments for the normalized 
bispectrum (8), under the assumption of Gaussianity. We assume without 
loss of generality l\ < 1% < 1$, and define 

f 1, for/i</ 2 <^3, 
A hl2l3 = 1 + 8\\ + 6ll + 3<$£ = I 2, for h = h < h or h<h = h, 

{ 6, for li = h = h; 

here and in the sequel, 5 h a denotes Kronecker's delta, that is, 6^ = 1 for a = b, 
zero otherwise. 

Under Gaussianity, it is obvious that the expectation of all odd powers 
of Ii x i 2 i z is zero. To analyze the behavior of even powers, we first recall 
that, for a multivariate zero-mean Gaussian vector (x\, . . . ,X2k)i we have 
the following diagram formula: 

(9) E(xi x x 2 x • • • x x 2 k) = ^2(Ex h x i2 ) x • • • x (Ex^^x^) , 

where the sum is over all the (2k)\/(k\2 k ) different ways of grouping (x±, . . . , 
X2k) into pairs (e.g., [1], page 108). Even powers of Ii^^ yield even powers of 
the aim's, which have a complex Gaussian distribution, weighted by Wigner's 
3j coefficients; we shall then need to use some arguments from graph theory, 
which is widely used in physics when handling Wigner's 3j coefficients (see 
[34], Chapter 11). 

Consider the Cartesian product / (8> J, where /, J are sets of positive 
integers of cardinality #(/) = P, #(</) = Q', it is convenient to visualize these 
elements in a P x Q matrix with P rows labeled by i and Q columns labeled 
by j. A diagram 7 is any partition of the P x Q elements into pairs like 
{(*l>ii)i (*2, J2)}; these pairs are called the edges of the diagram. For our 
purposes, it is enough to consider diagrams with an even number of rows 
P; we label T(I,J) the family of these diagrams. It can be checked that, 
for given I, J, there exist (P x Q — 1)!! different diagrams, each of them 

composed of (P x Q)/2 pairs; we recall that (2p - 1)!! = (2p - 1) x (2p - 

3) x • • • x 1 for p = 1,2, We also note that if we identify each row ik with 

a vertex (or node), and view these vertices as linked together by the edges 
{(*fe)ife)> (^ifc'i 3k')} = ikik', then it is possible to associate to each diagram 
a graph. As it is well known, a graph is an ordered pair {I,E), where I is 
non-empty (in our case the set of the rows of the diagram) , and E is a set of 
unordered pairs of vertices (in our case, the pairs of rows that are linked in a 
diagram). We consider only graphs which are not directed, that is, (iii?) and 
(z2«i) identify the same edge; however, we do allow for repetitions of edges 
(two rows may be linked twice), in which case the term multigraph is more 
appropriate. In general, a graph carries less information than a diagram (the 
information on the "columns," i.e., the second element jk, is neglected), but 
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it is much easier to represent pictorially. We shall use some results on graphs 
below; with a slight abuse of notation, we denote the graph 7 with the same 
letter as the corresponding diagram. 
We say that: 

(a) A diagram has a flat edge if there is at least one pair {(h,ji), (12, J2)} 
such that i\ = 12] we write 7 G Tp(I,J) for a diagram with at least a flat 
edge, and 7 G r-p(I, J) otherwise. A graph corresponding to a diagram with 
a flat edge includes an edge ij.ik which arrives in the same vertex where it 
started; for these circumstances the term pseudograph is preferred by some 
authors (e.g., [11]). 

(b) A diagram 7 G 1^(7, J) is connected if it is not possible to partition 
the i's into two sets A,B such that there are no edges with i\ G A and 12 G 
B. We write 7 G Fc(I,J) for connected diagrams, 7 G T-^t(I,J) otherwise. 
Obviously a diagram is connected if and only if the corresponding graph is 
connected, in the standard sense. 

(c) A diagram 7 G r-p(I, J) is paired if, considering any two sets of edges 
{(k,ji), (12, 32)} and {(13, J3), (k, J4)}, then h =i 3 implies i 2 = u; in words, 
the rows are completely coupled two by two. We write 7 G Tp(I, J) for paired 
diagrams. 

It is obvious that for P > 2 a paired diagram cannot be connected. Note 
that if Q is odd, paired diagrams cannot have flat edges, so that the as- 
sumption 7 G r-p(I, J) becomes redundant. 

(d) We shall say a diagram has a k-loop if there exists a sequence of fe 
edges 

{(k,ji), (i2,j2)}, ■ ■ ■ , {(ik,jk), {ik+i,jk+i)} = (^2), • • • , (ikik+i) 

such that i\ = ik+i] we write 7 G ^^^(1, J) for diagrams with a fc-loop and 
no loop of order smaller than k. 

Note that Tf(I,J) =r i ( 1 )(I, J) (a flat edge is a 1-loop); also, we write 

Tc m (i,J) = rc(i,J)nr m (i,j) 

for connected diagrams with fe-loops, and T c j^(I,J) for connected dia- 
grams with no loops of order k or smaller. For instance, a connected dia- 
gram belongs to T c jj^j(I,J) if there are neither flat edges nor two edges 

{(*l> 3l), (h, 32)} and {(13,^3), (u, 34,)} such that i x =i 3 and i 2 =u\ in words, 
there are no pairs of rows which are connected twice. 

A graph is Hamiltonian [written 7 G Th(I,J)] [11] if it has a spanning 
cycle, that is, if there exists a loop which covers all the vertices without 
touching any of them (other than the first) twice. Two graphs G\ = (I\,Ei) 
and G2 = (I2, E2) are isomorphic if there exists a one-to-one, onto mapping 
4>:I\ — > I2 such that i\i2 G E\ <^=^ 4>(ii)(j)(i2) G £2- 
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In many cases either I or J (or both) can be simply taken as the set of 
the first p or q natural numbers, that is, I = {1, . . . ,p}, J = {1, . . . , q}. Under 
such circumstances, when confusion is possible we shall occasionally write 
T(I,q),T(p, J) or T(p, q) for T(I, J). Some examples of graphs are drawn in 
Figures 1-6. 

We have the following result. 

Theorem 3.1. For all h<h< h we have 

(10) Elf ll2l3 =A hhl3 ; 
moreover, for p = 2,3, 4 ; 

(11) £?i* jB = (2p-l)!!A^+0(lT 1 ). 

Proof. Result (10) is known in the physics literature; see, for instance, 
[19]. For (11), we recall that 

(12) Ea v m . . oj. m . =(-l) m *i«Ci S l , j2 S~^T 2 ; 

hence, in view of (9), and because the spherical harmonic coefficients are 
(complex) Gaussian distributed, the following formula holds, for all /: 

as) 4nn%4= £ ^hJM 

U e /j=i V^h J 7G r(/,3) 

where we define 

(i4) 5( r ,h,h,h)= n (-i) miuj -'C!^^:- 



10 



D. MARINUCCI 



For any diagram 7, we can also write 

' h h h 

j=lrriij=—lj ) iG/ 



(15) %ii,uiHnn e ini; 4 ^j^^s), 



where 

Inn e }=e-e. oi>-,w=/ ; 

Ug/ j=lm,ij=-lj ) vn,i x \ m ip3 

to be more explicit, there are 3 x P summations to compute: for instance, 
when P = 2 we get 

f 3 h ~\ h h h h h h 

II II E = E E E E E E 

viG{ii,«2} i=lm t j=- 2j J ro» 1 i=— h mj 1 2=— i2 "i il3 =— ^3 mj 2 i=— £1 m i2 2=—h m i23 =—l3 

Furthermore, we also define 
D[T(I,3);h,l 2 ,h} 

= Y, D[ T) hMM] 
7er(/,3) 

={nn rJnU £ 

in words, !)[•; Zi, Z 2 , Z 3 ] represents the component of the expected value that 
corresponds to a particular set of diagrams. Notice that 



h 


h 


E • 


■ E 


mii=— h 


m,2 P ,3 = 


h 


h 


E • 


■ E 


m 1 i=-l 1 


"l2p,3 = 


h 




E • 


• E 


m 1 i=-l 1 


«l2p,3 = 



3 



l\ Z 2 I3 \ tt a ljn\ij 



2p i 

ma m i2 m i3 / ^ ^ 



' Zi Z2 Z3 
ma m i2 m i3 



) 7 6r(2p,3) 



D[r(2p,3);Zi,Z 2 ,Z a 



Now 



D[F(2 Pl 3);h,l 2 , h] = D[T P (2p, 3); h, Z 2 , Z 3 ] 

+ J D[r(2p,3)\rp(2p,3);Zi,Z 2 ,Z 3 ]; 
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our aim is to show that 



(16) 



D[r P (2p ) 3);Zi,Z 2 ,Z 3 ] = (2p-l)!!Af i 



hhh' 



(17) 



D[F(2p, 3)\r P (2p, 3); h, l 2 , l 3 ] = 0(1^). 



To establish (16) and (17) we rely on Propositions 3.1 and 3.2, whose proofs 
are collected in the Appendix. □ 

The next lemma relates to the "Gaussian" component of the expected 
value, that is, the diagrams that are paired. 

Proposition 3.1. For any p G N, and I with cardinality = 2p, we 
have 



The proof that (17) = 0{l l ) requires considerably more work; the next 
three lemmas refer to diagrams with loops of order 1, 2 and 3, respectively. 

Lemma 3.1. For diagrams with a flat edge, 7 G Tp{I,2>), we have 



The next two results show how diagrams belonging to Fc L r 2 \(I,3), 
T(7i(3)(/, 3) can be "reduced"; namely, they show how the corresponding 
summands in the expected value can be expressed in terms of smaller- 
order diagrams. Without loss of generality we can take #(/) > 4; indeed 
the case = 2 has been dealt with in the proof of Theorem 3.1, while we 
know that odd moments are identically equal to zero. Let 7 G ^cl(2){I^) 
be a connected diagram with a 2-loop, and denote by ii,i 2 the rows that 
are linked by two edges; in other words, the diagram includes both the 
edges [(ii,ji),(«2,Ji)] and Kh,32), («2, j' 2 )Y, m th e sequel, j k ,j' k takes values 
in (1,2,3), for any integer k. Because the diagram is connected, there must 
exist also edges [(n, j 3 ), (i 3 , j' 3 )\ and [(i 2 , j'4), («4, fd], where i 3 ,U7^i>*2- We 
denote by jR(i lt i 2 ) the lower-order diagram which is obtained by deleting 
[(*i,jl),(*2,Ji)] and [(ii,i2),(«2)J2)]) and substituting [(h, j 3 ), (i 3 ,f 3 )] and 
[(i 2 , ji), with [(i 3 ,j 3 ), {kJi)]. In graphical terms, 7R(u,i 2 ) is obtained 

by cutting the two nodes i\,i 2 and merging together the edges that departed 
from them to reach other vertices; 7ij(i 1) i 2 ) can itself belong to r^p) (I~ 2, 3) 
and the argument can be iterated (Figure 2). We note that these reductions 
need not be unique in general; any arbitrary choice of a suitable pair of 
nodes would not affect our argument, however. 



D[Tp(I,3);hJ 2 ,l 3 ] = (2p-iy.\Al 



D[ T ,h,i 2 ,i 3 ] = o. 
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Fig. 2. 7er CL(2) (6,3) and 7fl(n,* 2 )- 



Lemma 3.2. Forj£ ^cl(2){I^) an ^7_R(u,i 2 ) as defined before, we have 
D[r,h,h,h] = 2l +1 D[jR(i lii2 y,li,l2,h]- 

By definitions (14) and (15) D[y; l\, h, h] can be nonzero only if lj 3 = lj> = 
lj 4 = , so there is no notational ambiguity in Lemma 3.2. More explicitly, 
assuming, for instance, that the edges [(1, 1), (2, 1)] and [(1, 2), (2, 2)] are 
present in 7, then a factor (2/3 + l) -1 will emerge from the reduction. 

We now focus on diagrams with a 3-loop. Let 7 G ^cl(3) (^1 3) be a con- 
nected diagram with a 3-loop, and denote by i±, 12, 13 the rows that are linked 
by the loop; in other words, the diagram includes the three edges 

[(u,ji),(«2,j2)], [(ii,h),(h,h)\, 
Because the diagram is connected, there must exist also edges 

[(h, J», (u,js)], [{hih), (k,jio)}, [fe, J11), (kJn)] 

where ^4,25,^6 7^i^2^3- We denote by jRU 1} i 2t i 3 ) the lower-order diagram 
which is obtained by replacing 12,13 with i± and then deleting all flat edges. 
More explicitly, 7_R(i li i 2 ,«3) i s obtained by deleting 

[(n,ii),(«2,j2)], [{h,h),(h,h% [(*3, h), (h, je)], 
and substituting 

[(*2, Jq), (*5, Jio)], [(hJll), («6, Jl2)] 

with 

[(n,ig), («5) jio)], jii)) («6,j'i2)]- 
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Fig. 3. 7 g r CL(3) (6,3) and 7H(ii,t 3 ,i 3 )- 

In graphical terms, we are merging three nodes into a single one (see Figure 
3); again, jR(i lt i 2t i 3 ) can belong to T CL ^(I — 2,3) and the argument can be 
iterated. 

Lemma 3.3. For 7 G 1^(3) (I, 3) a connected diagram with a 3-loop and 
7/1(11,12,13) as defined before, we have 

D[T,h,h,k] = ^ ^ I ^[7/1(11,12,13)^1^2^3], 

where on the left-hand side we have used Wigner J s 6j coefficient, defined in 
the Appendix; hence 

D[r, h, h,k] = O^D^n^^y^h, h\). 

The proofs of Lemmas 3.1, 3.2 and 3.3 can be obtained as an application 
of the graphical method described, for instance, in ([34], Chapters 11 and 
12), and they are hence omitted for brevity's sake. 

The next proposition exploits the previous results to provide a bound 
(17) on the "non- Gaussian" part of the higher-order moments of the angular 
bispectrum. 

Proposition 3.2. For all I such that #(/) = 4, 6 or 8 and l\<h< k, 
we have 

D[T(I,3)\Tp(I,3y,l 1 ,l 2 ,l 3 }=0(l^). 
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Remark 3.1. A careful inspection of the proofs of Propositions 3.1 and 
3.2 reveals that for l\ < I2 < h we obtain the special case 

4 o , 6 6 6 , «Hi h h 



nR) EIhl ^~* + 2l 1 + l + 2l 2 + l + 2l 3 + l +6 \h h h 
^ lH > 6/1+9 

a result that we shall exploit for statistical applications. 

Remark 3.2. It is natural to conjecture that a result analogous to 
Proposition 3.2 will hold for all sets = 2p, p £ N. This conjecture is 
not simple to explore, however, because the proofs in the Appendix require 
some analytic properties of the summations of (products of) Wigner's 3j 
coefficients, which have not been extended (to the best of our knowledge) 
to products of arbitrary order. 

4. Unknown angular power spectrum. In this section we focus on the 
more realistic case where the angular power spectrum is unknown and es- 
timated from the data; so we consider rather than h^i^- As before, 
under Gaussianity of the underlying field T(9, ip) 

by a simple symmetry argument. Now note that 
\ a io\ 2 2\an\ 2 2\au\ 2 



Ci Q Q 

\a>io\ 2 2|ctnJ 



(2l + i; 



|a«)| 2 +Em=l 2 km| 2 Wiol 2 +Em=l 2 l 

2\au\ 2 



\ a lo\ 2 + Em=l 2 hr 

d ^ f (2Z + l)(6o,...,^) = (2Z + l)DirQ,l,...,l); 

here = denotes equality in distribution and Dir(0Q, . . . ,9 P ) a Dirichlet dis- 
tribution with parameters (9q, • • • 1 &p)- Define 

/,m a lm ~ a lm n 1 

(19) u lm = ^=, u im = —j=, m = 0,l, 

V w 



We have the following simple result. 
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Theorem 4.1. Let I and p be positive integers, and define 



g{hp) = tl 

k=l 



21 + 1 



2l + 2k-\\' 

Now for u and u defined by (19) , we have 



El UlQ... UlQ un... un u n ...u* x ...u lk ... U lk u* lk . . . u* k 



q times q 1 times ^ times 



q k times times 



E< uiq... uiq un... un u n ...u a ...ui k ...ui k u* k . . 



u 



Ik 



qo times gi times ^ t i mes 

x g(hqo + qi H h^). 



q k times q ' t i mes 



Proof. By symmetry arguments, it is easy to see that both sides are 
zero unless qo = 2po (say) is even and qi = q\=Pi (say), for i = 1, . . . ,k. The 
ui m are independent over different m's, and thus we have 



E 



uiq . ^ . uipUn • ■■ ■ unu n . . . u* lx ■ . . u ik . . ^.ui k u* lk . ..u* lk 

2p times gi times times Ik times 



<4 1 



Eu^Eu%(u h r...Eu%( 

k 

(2p -l)\ll[p i \, forp >0 
for po = 



i=i 



because 



Y[pi\g{l;p), 
(. i=i 



u w = N(0, 1) and u/ m «* m = |u; m | 2 = exp(l). 



Now write p = po + • • • + Pk, and note that ([18], page 233) 
(2/ + l)f 



RfPO . . . f:Pk 

2pi+-+Pk 

(2z + i)p r(/ + i/2) r(p + i/2)r(p 2 + i)---r( Pfc + i) 
2Pi+-+p*r(/ + p+i/2) r(i/2) 

{2p - l)\\p x \ x ■•■ xp k \ 



(21 + 1) x • ■ • x (2Z + 2p - 1 



-(2Z + r 
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i=i 

k 

~[pi\g(l;p), 



for po > 0, 
for po = 0, 



i=i 

as claimed. □ 

Some special cases are as follows: 
,ho| 2 



E 



Ci 

(\aio\ 2 " p 



C, 



m = +1, . . . , +1, 



(2p-l)!! 



E 



I "5 



(oi _i_ 

(2Z + 1) x ••• x (2l + 2p-iy ' ' 

(2/ + i)p r(z + i/2) r(i/2)r(p + i) 
2? r(z+p + i/2) r(i/2) 



p! 



(2Z + 1) x ••• x (2Z + 2p-l) 
and for p = pi +p 2 , pi,p2 > 



(21 + If, 



E 



^olVYM'Y 2 ! (21 + iy r(z + i/2) r( Pl + i/2)r(p 2 + 1) 



v di J \ di 



r(i/2) 

(2Z + 1)p 



an I 



2\ Pi 



2\ VI 



2P T(l + p + 1/2) 

(2p!-l)!!p 2 ! 
(2/ + 1) x • • • x (22 + 2p - 1) 

(2Z + 1)P r(Z + l/2) 
2P T(l + p+ 1/2) 

r(i/2)r(p 1 + i)r(p 2 + i) 
r(i/2) 

= (21 + 1)P 

(21 + 1) x ••• x (2l + 2p- l) v ; ' 

By Theorem 4.1, it is possible to establish a simple relationship between the 
normalized bispectrum with known or unknown angular power spectrum. 
More precisely, it is immediate to see that for l\ < Z 2 < h 



i=l 



that is, for instance, 
(20) 



EL 



EL 
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and 

(21) ^ fel8 = fai. (1 - ^ + 3) C 1 ~ 2fa + s) C 1 2/3 + 3, 

Also, for li = I2 < I3 and l\ < I2 = I3 

E %hk= EI U^9{h\2p)9{W,p), EI^ lzl2 =Elf^ h g(h;p)g(h-2py, 
finally, for l\ = I2 = I3 = I 

EI% = EI%g(l;3p), 

so that, for instance, 

(22) Eif ll = 6(l--^-](l 



2l + 3j\ 21 + 5 



It is interesting to note that 



EL 



(23) Ei% 2h <EI% 2h and Km-^ = 1 

for all choices of (h,h,h) f° r which the bispectrum is well defined. 

5. Some statistical applications. In this section we exploit the previous 
results to derive the asymptotic convergence of some functionals of the bis- 
pectrum array. Because the expected bispectrum is identically zero for Gaus- 
sian fields, these functionals arise as natural candidates in the development 
of statistical tests for non-Gaussianity. 

More precisely, assume that the resolution of the experiment is such that 
it yields a maximum observable multipole equal to L. It is, in practice, 
infeasible to take into account all available bispectrum ordinates for the 
implementation of a statistical procedure: indeed, for current experiments 
these ordinates are of the order of L 3 ~ 10 8 /10 9 , and the evaluation of all 
these statistics is beyond the power of the fastest supercomputers for the 
near future. It is therefore mandatory to consider only a subset of bispec- 
trum ordinates for the test. There are, of course, several possible choices of 
configurations. We shall restrict our attention to two of them; precisely, for 
finite integers Iq>2, K > we shall consider the processes 

1 [Lr] f 1 K h nil u 1 

(24) JiL-,i ,K(r) = —= { rrr—T X \ , 

(25) J2LI k{t) = 1 X] ( 1 X - 1 
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and 



1 [Lr]-l -K r 1 k \ 

(26) J3L;l ,K(r) = —/= { l T r i i E ^O+ti.M+Jo+u f' 

V I i =Zo+A '+i I V A + 1 u=0 J 

1 [i>]-/o--K' f -, K f2 

J 1 lo+u,l,l+l +u 



(27) j 4M (-) = 7f E TTfrrE 



where [■] denotes the integer part of a real number; < r < 1 and Iq is an 
(arbitrary but fixed) value which can be taken equal to 2, for instance, for 
cosmological applications (remember the dipole Iq = 1 is usually discarded 
from CMB data, as it is associated with kinematic effects mainly due to the 
motion of the Milky Way and the local group of galaxies). As usual, the 
sums are taken to be equal to zero when the index set is empty. K is a fixed 
pooling parameter: for K = we obtain the special cases 



i [Lr] f 



(28) 



and 



(29) 



JlL;l {r) = -= J2 

V L / 2 lcven , >lo V6 



[Lr 



1 [Lr]-l ^ 

hm {r) = —i= 2J J io,i,i+io- 
VL «=«o+l 

[Lr]-i f2 

1 X-^ 1 lr 



=io- 

The normalizing factors are chosen to ensure an asymptotic unit variance 
for all summands, by means of Theorem 3.1 and (18), (21) and (22). For 
instance, 

Vax{ if 0tl>l+l0 - 1} 

= Var{Ij t ij + i } = EI lQjU+lo — {EI loll+lo } 
4 / 2Zp + l \ / 2H-1 \ / 2Z + 2M-1 

( 30 ) khl ' l+lo \2l + 3j \2l + 3j \2l + 2/o + 3 

2/ + l 

,2/o + l * ' l2/o + 3 

2 + Oil- 1 ). 



t + 0/ \zl -t- zto + >jy 



The two pairs of processes JiL;i ,K(r), J 2 L-,i M r ) and J 3L;l ,K(r), ^4L ; i ,^( r ) 
can be viewed as sorts of boundary cases for the possible configurations 
of multipoles. Indeed, although none of them has so far been considered 
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in the literature (to the best of our knowledge), it seems natural to view 
JiL-,i ,K( r )i J2L;i ,K(f) as proposals very close to much of what has been done 
so far in CMB data analysis. More precisely, it has become very common 
to restrict attention to multipoles close to or on the "main diagonal" l\ = 
li = Is = I, under the unproved conjecture that the greatest part of the 
non-Gaussian signal should concentrate in that area. This same rationale 
motivates Jxl\1 ,k{t), J2L;1 ,k(t)\ we shall show below, though, how such a 
choice can be very far from optimal in relevant cases. On the other hand, 
J3L;1 ,k(i~), ^4L;Z ,-^( r ) rer y 011 a sor t of opposite strategy, that is, for a fixed Iq 
we aim at maximizing the distance among multipoles, albeit preserving the 
triangle conditions li < L + Ik ■ There are several alternative procedures one 
may wish to consider, but those we mentioned lend themselves to a simple 
analysis, while highlighting some quite unexpected features of asymptotics 
for fixed-radius fields. 

Theorem 5.1. As L — > oo, for any fixed integers Iq > 0, K > 0, 

JlL;l ,K(r), J2L;l ,K(r),J3L;l ,K{r), JiL;l Q ,K{ r ) W(r), < r < 1, 

(31) 

where denotes weak convergence in the Skorohod space D[0, 1] and W(r) 
denotes standard Brownian motion. 



Proof. The proofs for the processes J a L\i .K{ r )-, a = 1, ... ,4, are very 
similar; we give the details only for the most difficult case, namely J4,l-,i ,k{t)- 
Here the proof is made harder by the complicated structure of dependence; 
note indeed that the set of random coefficients {ai m :l = Iq, . . . ,Iq + K,m = 
belongs to each summand in (29). Denote by 3 s ; the filtration 
generated by the triangular array {ai-i, . . . , aij}, I = 1,2, ... , and define 

X t , L = - 7 ^=f; ^ +,, ^ + "" 1 > l = l + K + l,l + K + 2,..., 
v K + 1 u =o v 2L 

that is, 

[Lr]-l ~K 

■hL;l ,K(r) = X hL- 
l=l +K+l 

Now we note first that 

(32) = w+th ^ 

= 0, m>l, 
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because for all < l\ < I < I2 < I3, 



2 / h h h \ ( h h h 



E{if lhl3 m= E E 

7711,7712,7713 rn' 1 ,m' 2 ,m' 3 



mi m2 1713 J V m l m 2 m 3 



x E 



a hmi a l 1 m' 1 ^limi a lim' 2 a hrn,3 a l 3 m' 3 



ChCl 2 Ci 



E E 



mi, m 2 , m 3 m' 1 ,m' 2 ,m' 3 



h h h \ ( h h h 

mi 1712 m 3 J \ m l m 2 m 3 



a lirn 1 0'l 1 m'i ,_, f a l2rn2 a l2rn' a l 3 m 3 a Um' 
X -E> - 



Ci I Ci 2 Ci 3 
h h h\( h h h \ a 'i'»i a 'iS 



E E 



1 \n |2 

E 1 \ a hmi\ _ * 
2h + l Ci 

mi 1 

Equation (32) does not imply that the triangular array {-2Q,l}/=2,3,... obeys 
a martingale difference property, because the sequence X\^ is not adapted to 
the filtration 3?/. However, (32) proves that the pair sequences {X^l, ^1)1=2,3, .. 
do satisfy a mixingale property [7, 25], that is, 

(33) [E{E{X^ L \%_ m }) 2 \ 1 / 2 < ci^-t for m> 1, 

V L 

(34) [E(X ltL - E{X hL \Qi +m }) 2 ]^ 2 < C2^-jt for m > 1, 

(35) for some ci, C2, </>> 0. 

Actually the left-hand sides of (33) and (34) are identically zero for m larger 
than Iq + K, so that, for suitable choices of the constants ci,C2, the bounds 
on the right-hand sides hold for an arbitrarily large </>. Note that {X? L ,l = 

1, 2, . . . , L, L = 1, 2, . . .} is a uniformly integrable set, because if { i 2 i 3 has finite 
fourth-order moments which are uniformly bounded [Theorem 3.1 and (23)]. 
Also, it is readily seen that 

sup iimsup ■ — ; < co, urn max EXi L = 0. 

0<ri<r 2 <l L^oo ^2 ~ 1 L-+00 i=l,...,L 

We have thus established conditions (2.2), (2.3) and (2.5) in [25] for the 
functional central limit theorem to hold. To complete the proof, we only 
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need to show that 



lim E 

L—KX> 



( / [Lt] \ 

E \ E Jfw 

l \l=\Ls] J 



\hr\ \-(t-a) 



for any r < s <t. 



For notational simplicity and without loss of generality, we consider the 
special case K = 0. We have 



E< 



/ [it] v 

Vz=[Lsl / 



3[Lr] 



[Lt] [Lt] [Lt] 

£ £{Xj L |3 [£r] } + 2 X) E ^{^,L^,L|9[Lr]}- 
J=[La] l=[Ls]l'=l+l 



Now for 1,1' > [Lr] we have 

-EH^L^Z'.Ll^Lr]} 

1 



WT E {( I l ,l,l+l ~ 1 )(^ , 



(36) 



2L 
1 



l',l'+l 



1)IW 



~ E{I? 0!lil+lo \% [L r]} ~ #{i? oW d%r]} + 1] 

2L 



There are now two possible cases, namely I + Iq = V and I + Iq ^ V; in the 
latter 



(37) 



E { I lo,l,l+l I lo,l',l'+lo\^lLr]} 

E E 

™0l ,«i02 ,mo3 ,"i 4 mu ,mi2 ,m 2 i ,m 2 2 

x Mo I l + k 

\m 02 -mn -m 2 i 

x Mo Z' + Z 

Vm 03 mi2 m 2 2 . 

a lomoi a lomo2 a lomo3 a lo'mo4 



Mo i M/ 
v m i mn m 2 i 



Zo I' I + Zo 

m 4 —77112 -"l 2 2 



Cfo 



E 



moi,mo2,"io3,mo4 



^m i° 2 ^-mo 3 04 a lom i a l m 2 fl tom 3 a lom 04 

(2Z + 1) 2 ' C» 
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the second to last step following from (49) in the Appendix. Otherwise, for 
l + la = l', 



(38) 



E E 

moi ,m 02 ,m 03 ,m 4 mil ,mi2 ,m 2 i ,m 22 
/ l I 

\m 02 -m u 
( k I' 



X 
X 



( h I I' 
\m i mn m 2 i 

Iq I' I' + Iq 



m 3 m w m 2 2 



-77121 

^ + ^0 i a lomoi a lomo2 a lomo3 a lo'mo4, 

-m 2 2, 



+ 2 



E 



E 



m 01 ,m 2 ,m 3 ,mo4 mil ,mi 2 ,m 2i ,m 2 2 

Z I I' 



Vm i mu m 2 i 

/ / V l' + h 
m 2 m u -m 2 i J \m 03 -mn m 22 

Iq I I + k \ a lom iO'l m 2 a lom 03 a l m 4 



m 4 -mu -m 22 



1 + 2B, 



where 
5 = 



E E 

moi,mo2,m 3,mo4 mii,mi2,m 2 i,m22 
m 03 -mn m 22 J \m 04 



Iq I V 
m i mu m 2 i 



l I 

m 02 mu 



I' 
-m 2 i 



x 



V I O-l m O i a lom 2 a lom 03 ( ll mo4 

-mi2 -m 22 . 



r r 2 
°z 



s=0a 
X 



Iq I + Iq I' 
l s Iq 



E 



( k 



m i,m 2,m 3,m 04 



( k 



\m i m 02 crj \m 03 m 04 a 

a l m 01 a l m 2 a lom 3 fl lomQ4 
x ^ , 

in view of (54) in the Appendix. Now the triangle conditions entail that the 
summands are nonzero only for s < 21q; from (50) and (55) below we learn 
that 



k k s 
, m i m 02 a 

Iq I I' 

k s l 

whereas it is also trivial to note that 

|2 



f k k s 

V m 03 JTi 04 cr 

) l' + k I' 

) s Iq 



< 



1 



< 



y/2l + 1 ' 
1 



< 



1 



a lomo 1 a lomo2 



lo 



< max 



\ a lor 



a 



(2Zq + 1) max 



V2U + 1 ~ V / 2T+T' 

\n I 2 
| u Zom| 



lo 



J2k=-l \ a lok\ 



< 2l Q + 1. 
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(39) 



l*l<E E 



(2Z + 1) 



E 



1 



< 



9/4-1 Z ^ 9/ I 1 

s<2l \<t\<2s m i,m 2,mo3 3 mo4 

(4Z + 1)(2Zq + 1) 7 
2Z + 1 



x (2Z + if 



Combining (36), (37), (38) and (39), we have 



\E{X l>L X VtL \% [Lr] }\ < 



1 (4l + l)(2l + l) 7 



L 



21 + 1 



J l+lo> 



[Lt] [Lt] 

E E{X ltL X VtL \Qt [Lr] } 

l=[Ls]l'=l+l 



< 



logL(4Z + l)(2Zo + l) 



L 



21 + 1 



The remaining part of the proof is similar. More precisely, we have 
E{X? L \Z [Lr] } 



-Z E {Il ,l,l+l ~ M^[Lr]}, 



E {Il ,l,l+lo\^[Lr}} 

= 3 E 



E 



moi .moa ,i«03 ■"'■04 mn,mi2 ,m2i,m 2 2 











Mo i Wo 
V?n i mil "i2i 



Mo ^ Z + Zo 
\m 2 -mn -m 2 i 

/ l + l 



mo3 mi2 m 22 / \mo4 -m i2 -m 22 

a lomoi a lomo2 a hmo3 a Zomo4 



+ 6 



E 



E 



"H)l,«t()2. "'03, "104 mil ,mi2,"l21,?"22 







I l + l 



Iq I I + Iq 

m i mn m 2 i 

Mo ^ l + h 

\m 02 -mn m 22 

/ l + h 



m 3 m i2 -m 2 i/ \ m 04 -mi 2 -m 22 

a £o"K)i a hmo2 a lomo3 a lomo4 



cf 



3A + 6B. 
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As before, 



E 

raoi,mo2,mo3,mo4 

1 l„. 



(2Z + 1)2 " 



f' 2 



E 



| a Zo"iOll l a £omo3l 



moi,mo3 



(2/ + 1) 



and [see (54)] 



E 

mii,mi2,rrt2i,m22 
Z / 



?o i l + k\ f h I l + k 
m i m,n m 2 i / \ m 02 -"in rn 22 



/ 



m os mi 2 -mix) \ m 04 -mi 2 

E £<*•+!>{£ 'V 1?}' 



S<2Zq cr=-s 



^0 



m i m 02 cry \ m 03 mo4 <7 



Z + Z 
-m 22 



o l o 



whence 



B< 



(4Z + 1)(2Z + 1) 7 



21 + 1 

as in (39). We have thus shown that, for some C > 0, 

[Lt] X 2 



/ +t\ \ ■ 

E 

\l=\Ls] 1 



\lt\ > - (t - a) 



< 



[Li] 

E (£{^ 2 L i% r] }-i; 

J=[Ls] 



+ 2 



[Lt] [Lt] 

E E £{^,L^',L !%,-]} 
l=\Ls] l'=l+l 



+ 



L 



< 



C 



[Lt] 

E 

l=[Ls] 



(4Z + 1) 



2Z + 1 



+ logL + l 



<C(4/ + l) 8 ^-0 

Li 



as L — > oo. 



Thus the proof is complete. □ 



Theorem 5.1 can be immediately applied to derive the asymptotic dis- 
tribution under the null of several non-Gaussianity tests. For instance, we 
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might focus on sup 0<r<1 J a L;i ,K(r); by the continuous mapping theorem we 
obtain 

lim P\ sup J a L;l ,K{r) < x \ = P\ sup W{r)<x\ 

(4U) L^oo Lo<r<l J lo<r<l J 

= 2$0)-l, x>0,o=l > 2,3,4 ) 

where <&(•) denotes the cumulative distribution function of a standard Gaus- 
sian variable (for the last equality, see, e.g., [5]). We shall now discuss the 
behavior of these procedures under some examples of non-Gaussian spher- 
ical fields. The behavior of higher-order angular power spectra under non- 
Gaussian alternatives is an extremely important research topic in modern 
cosmology, and still almost completely open for mathematical research. Very 
few analytic results are available, whereas the cosmological debate is still 
open on the nature of the non-Gaussianity to be expected. A simple and 
popular model for non-Gaussian temperature fluctuations reads 

(41) T NG (0, ip) = T(9, if) + /nl{T 2 (0, <p) - ET 2 (9, V )}- 

as before, we take T{Q,ip) to be an isotropic Gaussian field with zero mean. 
For "small" values of the nonlinearity parameter /nl, (41) can be viewed 
as a general approximation for random fields with minor departures from 
Gaussianity: the quadratic term can be regarded as the leading factor in a 
Taylor expansion of a general field g(T(9, tp)), for a suitably regular function 
g(-). Equivalently, the terms on the left-hand side can be considered to be 
the first two elements of expansion of g(T(0,(p)) into a series of orthogonal 
Hermite polynomials H q (T), q = 1,2,.... For these reasons, (41) is very 
widely adopted to represent the primordial field of temperature fluctuations 
in cosmological models of inflations, which stand now as the leading models 
for the dynamics in the primordial epochs around the Big Bang [3, 19, 28]. 
It is known in the physics literature [19] that the bispectrum of (41) can be 
approximated by 

(42) B h i 2 i :i = G/NiAiZaZs ( q q q J {ChQ 2 + Ci 2 Q a + C h Ci 3 }, 
where G is a positive constant, 

/(2/ 1 + l)(2/ 2 + l)(2/ 3 + l)\ 1/2 

K l2h - y ^ ) , 

and lower-order terms are neglected. Expression (42) is known as the Sachs- 
Wolfe bispectrum; we shall take (42) as a benchmark model for our discus- 
sion of non-Gaussianity. We introduce a very mild regularity condition on 
the angular power spectrum, that is, we shall assume that C/ is such that, 
for fixed Iq> 0, 



(43) Ci+ioOcC,, 
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where oc denotes that the ratio of the right- and left-hand sides tends to 
a positive constant, as I — > oo. The assumption described by (43) is not 
unreasonable: indeed, C\ oc l~ a (for some positive constant a) for many, if 
not most, cosmological models. 

For simplicity, let us assume that the normalizing angular power spectrum 
is nonrandom, that is, known a priori; without loss of generality we take 
K = 0. Recall that ([34], equations 8.1.2.12 and 8.5.2.32) 



h h h 




(_l)(Ji+fe+J3)/2[( Zl+l2 + l3 ) /2 ]! 



Thus, for fixed Iq > 2, 

lo I I + Iq 




[(h +l 2 - h)/2]\[(h -h + h)/m(-h + h + h)/2]\ 

{h+h-hnh-h+hn-h+h+h) \ \ 1/2 

(-l)'o+ ; (/ + / ) X ••• X (l + l) 

w 



y/(2l + l) x ••• x (2Z + 2Z + 1) 



C 



Vi 



for some C > which depends on lo but not on I. Then we have easily that 

[Lr] 



(2Z + 1)(2Z + 1)(2Z + 2Z + 1) 



io+l 



4tt 



(44) 



Vl 



Ci Ci . Ci Ci+i CiCi + i Q 



C, 



+ 



[Lr 



OC 



JNL n CiCi /nl n , T 

V L i=i 0+ i V c Wo VL 1=1q+1 



+ 

[Lr] 



Likewise 



[Lr] [Lr] 

£;j 4L ( r)K J2 {Elt l , l+lo -l}>^=J2{( EI W+io) 2 -^ 

Of course, both (44) and (45) diverge as the number of observed multipoles 
increases (L — ► oo), that is, as the resolution of the experiment improves. 
The constants of proportionality are typically small; for the model we adopt 
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in the simulations below, they are of the order 10 4 for (44) and 10 8 for 
(45). 

On the other hand, for li=l2 = h = l, 

I I 1)S i/2.[3//2]! / [If \W 



0J v ' [(Z/2)!] 3 l(3Z + l 



r 3 fl(3I+l)/2 f ;3/+3/2 i 1/2 /_j_\3«/2 



oc (-1)31/2 L^r _± ; 

°^ L > /3/2(/+l) \ (3/)3(I+l/2) / 



OC 



where we have used Stirling's formula n!/(\/27rn n+1//2 e n ) = 1 + 0(Yln 1 ). 
Now recall that 

which implies IC\ = o(/ _1 ) (in realistic experiments, a weighting factor Gi 
should be included to account for the beam pattern of the antenna). Hence 
we have 



EJlL;l ,K{r) 

(46) 




/ l\ Cf 



0, 
as L — > oo. 



A similar heuristic argument can be used for J2L-,i Q ,K( r ), suggesting that 
these statistics may have very little, if any, power against alternatives of the 
form (42). 

To validate the previous heuristics, we present a small Monte Carlo study 
on the power of these testing procedures. To this aim, we generated 200 
spherical Gaussian fields according to the currently favored scenario for CMB 
fluctuations, the so-called ACDM model; we omit the details, which can be 
found, for instance, in [6]. Nonlinearities were then introduced by means 
of (41), which represents the simplest non-Gaussian model, as argued ear- 
lier. To ease comparisons with existing procedures, we follow the standard 
parametrization used in the astrophysical literature, yielding a variance for T 
of the order of Var(T) = ET 2 ~ 10~ 8 . We can then present a rough relation- 
ship between the value of the nonlinearity parameter /nl and the relative 
amount of the non-Gaussian signal, namely 

x/Varl/NL^} 



(47) v v ;; u ;^ J = * /NL x 10- 

v/VarjT} v 

We consider /nl = 0, 100, 300, 1000 and we focus on the statistics 
sup <r<i^3L;i ,^( r )> su Po<r<i ^4L;i ,^( r ) for ^0 = 2 and if = 0,2,4; see (40). 
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We omit reporting results for sup < r < 1 Jil;1 ,k{t), sup^^ J 2 L i i ,K(r) be- 
cause the power related to these procedures turned out to be negligible, 
as expected from (46). Other types of statistics may be considered with- 
out altering the main conclusions we are going to draw here. We start from 
the empirical sizes (type I errors), which are reported in Table 1; we write 
S a ,L for sup 0<r<1 J a ,L:2,K{f)-, a = 3, 4. We take L = 250, 500; these values are 
conservative: we recall that L is reckoned to be of the order of 600/800 for 
WMAP and 2000/2500 for Planck. Note that all values in Table 1 and those 
in Tables 3-5 are expressed as percentages (%). 

Results in Table 1 suggest that the asymptotic theory presented in The- 
orem 5.1 provides a good approximation for the finite-sample behavior in 
the case of J3,l; the approximation is slightly less satisfactory for J^l, but 
the results improve markedly with L, which is reassuring. Because the test 
statistics are free of nuisance parameters, it is also possible to derive directly 
threshold values under the null of Gaussianity by Monte Carlo replications. 
We adopt both approaches to derive the power properties reported in Tables 
3-5; the Monte Carlo critical values are reported in Table 2. 

Note how the tabulated values approach the asymptotic results [i.e., 1.645, 
1.96 — see (40)] when L increases. To analyze the power of the test, we 
consider /nl = 100, 300, 1000; from (47) we can argue heuristically that these 
values correspond approximately to 1%, 3% and 10% of non-Gaussianity in 
the maps, respectively. In Tables 3-5, T and A denote the power with respect 
to tabulated and asymptotic critical values, respectively. 

At first sight, the results reported seem quite encouraging, as compared, 
for instance, to existing methods such as the empirical process, wavelets 
and local curvature: see, for instance, [6] for numerical simulations on the 
performance of these procedures. We stress once again, however, that this 
comparison could be to some extent misleading, insofar as in this paper we 
are sticking to some simplifying assumptions that are unrealistic for CMB 
experiments (namely, the absence of gaps and noise in the observed maps). 
The results we report, however, certainly suggest that further investigation 
of these procedures under more realistic circumstances is worth pursuing. 
We also note that the statistics based on J3L;i ,K( r ) substantially outperform 
those based on J4,L-i ,K{f) f° r this range of parameter values; indeed, for the 
latter the power is nonnegligible only when /nl = 1000, that is, when the 
level of non-Gaussianity in the map is approximately 10%. However, a com- 
parison of (44) and (45) suggests that the power for J4L;i ,K(f) may improve 
more rapidly as the resolution of the experiments grows; this is an impor- 
tant factor to keep in mind, as it is expected that the satellite Planck will 
achieve observations with L of the order of 2500. Also, because J4,l-,i ,k{t) 
does not depend on the sign of the sample bispectrum, it can be expected 
to be more robust against other types of non-Gaussian behavior. It is also 
worth noting that the power of the tests grows with the pooling parameter 
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Table 1 

Empirical sizes (/nl =0), a = 10% (5%) 





K = 


K = 2 


K =4 


S3, 250 


9.5 (4.5) 


9.0 (4.0) 


8.0 (2.5) 


S3, 500 


10.5 (4.5) 


9.5 (4.5) 


9.5 (4.0) 


£4,250 


5.5 (3.5) 


3.0 (1.5) 


1.0 (0.5) 


£4,500 


6.5 (4.5) 


8.0 (4.0) 


6.0 (2.0) 



Table 2 

Monte Carlo critical values, 01 — 10% (5%) 





K = 


K = 2 


K =4 


£3,250 


1.61 (1.81) 


1.55 (1.83) 


1.61 (1.72) 


£3,500 


1.69 (1.90) 


1.63 (1.92) 


1.61 (1.80) 


£4,250 


1.44 (1.71) 


1.27 (1.51) 


1.06 (1.35) 


£4,500 


1.56 (1.88) 


1.52 (1.85) 


1.38 (1.67) 



Table 3 

Rejection rates for /nl = 100, a — 10% (5%) 





K = (T) 


K = 2 (T) 


K = 4 (T) 


K = (A) 


K = 2 (A) 


K = 4 (A) 


£3,250 
£3,500 
£4,250 
£4,500 


18.0 (10.5) 
20.5 (13.5) 
9.0 (5.0) 
10.0 (5.0) 


20.0 (12.0) 
30.5 (20.0) 
10.5 (4.5) 
10.0 (6.0) 


18.5 (17.5) 
39.5 (33.0) 
10.0 (4.5) 
10.0 (6.0) 


17.5 (7.0) 
20.5 (12.5) 
5.5 (3.0) 
7.0 (4.0) 


16.5 (9.5) 
30.5 (20.0) 
3.0 (1.5) 
7.0 (4.0) 


17.5 (11.0) 
38.5 (26.0) 
1.0 (1.0) 
6.0 (2.0) 






Table 4 
Rejection rates for /nl = 30C 


', a = 10% (5%) 






K = (T) 


K = 2 (T) 


K = 4 (T) 


K = (A) 


K = 2 (A) 


K = 4 (A) 


£3,250 
£3,500 
£4,250 
£4,500 


31.0 (26.0) 
47.5 (43.0) 
10.5 (5.0) 
12.0 (6.0) 


47.5 (38.5) 
88.5 (69.0) 
11.5 (5.0) 
11.5 (7.5) 


56.5 (51.5) 
88.5 (86.5) 
9.0 (5.5) 
11.5 (7.0) 


31.0 (21.5) 
48.5 (40) 
6.0 (3.0) 
9.0 (4.0) 


43.0 (33.0) 
88.5 (68.5) 
3.5 (1.5) 
9.5 (6.0) 


56.0 (41.5) 
88.5 (79.0) 
2.0 (0.5) 
7.5 (2.5) 


Table 5 

Rejection rates for / NL = 1000, a = 10% (5%) 




K = (T) 


K = 2 (T) 


K = 4 (T) 


K = (A) 


K = 2 (A) 


K = 4 (A) 


£3,250 
£3,500 
£4,250 
£4,500 


80.0 (73.0) 
98.0 (98.0) 
15.5 (10.5) 
37.5 (29.0) 


99.5 (99.5) 
100 (100) 
26.0 (18.5) 
57.5 (50.5) 


100 (100) 
100 (100) 
32.5 (18.5) 
71.5 (58.5) 


80.0 (70.5) 
98.0 (97.5) 
12.0 (7.5) 
34.0 (24.0) 


99.5 (99.5) 
100 (100) 
10.5 (5.5) 
54.0 (47.0) 


100 (99.5) 
100 (100) 

8.5 (3.5) 
60.5 (47.5) 
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K; for values larger than 6-8, however, the effect is nearly negligible and 
the computational cost becomes prohibitive. 

In our opinion, there are two remarkable features that emerge from this 
section. Equations (44) and (45) suggest that a testing procedure in har- 
monic space can yield consistent tests of Gaussianity even for fixed-radius, 
nonergodic fields, at least under the simplifying assumption of this paper. 
This is to some extent an unexpected result. The second remarkable fact is 
the huge impact of the choice of combined angular scales on the expected 
power under non-Gaussian alternatives. It is noteworthy that the common 
choice of a (close to) "main diagonal" configuration can yield negligible 
power, the expected value of the non-Gaussian signal decreasing to zero as 
the resolution of the experiment improves. The determination of the triples 
of angular scales (h,h,h) where the largest part of the non-Gaussian signal 
is to be expected, for a given class of models, represents an issue of great 
importance for future cosmological data analysis. 

6. Comments and conclusion. It is important to make clear that the 
asymptotic theory presented in this paper is of a rather different nature with 
respect to what is usually undertaken for random processes or fields. More 
precisely, we are not assuming that the information grows in the sense that a 
larger interval or region of observations becomes available, but rather we as- 
sume that the same region (spherical surface, in our case) is observed with 
greater and greater resolution. We labeled this framework high-resolution 
asymptotics, whereas the more standard case where the observed volume 
grows can be termed, as usual, large-sample asymptotics. The idea that 
some consistent inference can be drawn from a process observed on a fixed 
finite region is certainly not new; see, for instance, [32] for a fixed-domain 
asymptotics approach to study optimal linear prediction of spatial processes 
(kriging). The term infill asymptotics is also used in the statistical litera- 
ture with the same meaning. We believe that this paradigm will become 
more and more fruitful in the years to come, with many possible contexts of 
applications. In the case of cosmological research, a proper understanding 
of the nature of the asymptotic theory involved is likely to be quite rele- 
vant from both the theoretical and the practical point of view. In fact, note 
that sequential experiments to measure CMB radiation result in exactly 
the same last scattering surface being measured, while the resolution im- 
proves steadily over time. For instance, the above mentioned NASA satellite 
WMAP, launched in 2001, is observing the same surface as the ESA mission 
Planck, due to be launched in 2007: in terms of the standard large-sample 
asymptotics, no improvement should be expected. On the other hand, for 
statistical properties that can be consistently investigated as the resolution 
of the experiment grows, Planck does offer substantial new information, its 
expected resolution outperforming WMAP by a factor 3 or 4 (in any case, 
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Planck will offer other improvements besides better angular resolution, e.g., 
better polarization measurements). It seems therefore quite relevant to sug- 
gest, as we did in this paper, that Gaussianity tests may exist that are high- 
resolution consistent, at least under some simplifying assumptions. The fact 
that consistent inferences can be drawn from nonergodic random fields de- 
fined on a bounded domain has other important consequences if we focus 
on epistemological issues. The status of cosmology as a science is occasion- 
ally questioned, on the grounds that it is, in a sense, a discipline based by 
definition on a single observation (our Universe). The possibility to draw 
consistent inferences for fixed-radius random fields provides, in our view, 
a strong argument to consider the corresponding physical properties fully 
within the domain of scientific investigation. It is a challenging task to char- 
acterize, under general conditions, the complete set of properties on which 
high-resolution consistent inferences can be drawn. 



A.l. The properties of Wigner's coefficients. In this paper, we make 
extensive use of Wigner's 3j coefficients, which are a very powerful tool 
to represent properties of random fields which are invariant to rotations. 
These coefficients were introduced in the framework of the quantum theory 
of angular momenta; they are also widely used in algebra in the framework of 
representation theory. In this section we shall recall some of their properties 
for convenience. 

Wigner's coefficients are defined implicitly by 



Many explicit representations are also available, but for general values of 
li,nii they are lengthy and hardly informative. For instance, it can be shown 
that ([34], expression 8.2.1.5) 



APPENDIX 




Yi imi {9, tp)Yi 2Tn2 (9, ip)Y hm3 (9, ip) sin 9 dip d9 





(-1) 



x 




(h + mi)!(Ii - mi)!(I 2 + m 2 )!(/2 - m 2 )! 




(-l) 2 (/ 2 + h + mi - z)\{h - mi + z)\ 



z\(l 2 + h-h~ z)Kh + m 3 - z)\(h -h-m 3 + z)\ ' 
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where the summation runs over all z's such that the factorials are nonneg- 
ative. We list here some important properties, and refer to [34] for proofs 
and further discussion: 

(a) Wigner's 3j coefficients are real valued; 

(b) they are different from zero only if mi + ni2 + 7773 = 0; 

(c) (parity) for any triple li,fa, h 

h h h 



( h h l 3 \ = (_i)h+h+h ( 
\m\ iri2 7773 J V 



-mi — 7712 — 7713 



(d) (symmetry) for any triple li,fo,h 

h h h\ = (h k h\( k h h 
m\ rri2 TO3 / \ m 2 m 3 m i ) V^3 m i "12 

= (_i\ii+h+h ( h h h 

\77l3 771-2 Wl 



("I) 



u 2 +h ( h h h 
\ mi m 3 



fWi+h+h ( h h k \ 

\ 777,2 777-1 "13 / 



(e) (orthonormality) for any triple h,h,h 
h h h 

E E E 

m 1 =-l 1 m 2 =—l2 m 3 =—l 3 

and 



( h h h 

\ mi 777-2 m 3 



(h h L\(h l 2 L'\ 



(49) E E 

mi=- hm2=-h 
(f) (upper bound) for any h,l2,h 



\ mi m 2 M J V 7771 Tl2 M ' / 2L + 1 ' 



f * 2 * 3 ) =0([max{/i,Z 2 ,W]" 1/2 ) 

\ 777i 7772 7773 / 

(gl) (sum of coefficients, I) for any positive integers a, b 

(51) EC- 1 )"* (« - a ^)=(-l) a v^+T^; 

(g2) (sums of coefficients, II) for any positive integers a,b,c,d,e and / 

£ £(-ir +/+c+ *(° J :)(; J _ e , 

(52) x( a d f \( c b f 



a 6 —(f> J \7 /3 (/> 
a b e 
c d f 
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(g3) (sums of coefficients, III) for any positive integers a, b, c, d, e and / 
(see [34], 8.7.3.12) 

X^(--\\$+"/f a b c \ ( a f d \ ( e b d 

^ { ' \-a -(5 jj\a <p -6J\-£ (3 -S 
(53) «.ftJ v \ , 

- (--\\ b -t ( c f e \\a b c 

1 ' \l <P -e)\e f d 

(g4) (sums of coefficients, IV) for any positive integers a,b,c,d,e and / 
(see [34], 8.7.4.20) 

a b c Wd f e \fg b e\f j f c 
a -f3 -<y) \6 -cp -e ) \r) (3 e) \n tp 7 



E 

(54) = f_l\a-b+c+d+e-f 

a s j \ ( 9 s d \ ( b c a\ ( b e g 



x 



E 



-SO- 



CK a —fi J \7) a —5J\j s f J \d s f 



Equation (52) can be used as the definition of Wigner's 6j coefficient, 
which appears on the right-hand side (with curly brackets). We refer again 
to ([34], Chapter 9) for some (extremely complicated) explicit expression 
for these coefficients and many of their properties. For our purposes it is 
sufficient to recall the following: 

(h) for any positive integers a,b,c,d,e and / 



(55) 





b 




< min^ 






;} 


e 





1 



V(2c+1)(2/ + 1)' V(2o + l)(2d+l) 

I 



V(26 + l)(2e + l)/' 
(i) the 6j coefficient is invariant under any permutation of its columns. 



A.2. Proofs of technical lemmas. 



Proof of Proposition 3.1. Let {Ii,...,I p } be a partition of / into 
2x3 matrices made up with two of its rows; write T for the class of these 
partitions, which has cardinality (2p — 1)!! (the number of combinations by 
which we can match 2p rows two by two to form p pairs). It suffices to notice 
that 

v 

D[Tp(I,3);h,l2,h]= J2 HDppilijyjMs]- 

{h,...,i p }eii=i 
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For any ij, D[Tp(Ii, 3); l\, l 2 , h] is easily seen to include A/ 1 j 2 2 3 nonzero sum- 
mands, each of them of the same form, up to a relabeling of the indexes; 
more precisely, 

D[Tp{I h 2,);hMM 

= A hhh Y (-l)mn+m i2 +m i3 (h h h\ 

123 V ran m i2 rrti-x J 

mn,mi3,mi3 v 11 lA ™ ' 

h h h 
-ma -m i2 -m i3 

1 i 6 \ m m m 1 

m a ,mi3,m i3 v " lz la ' 

Thus 

D[T P (I,3);h,l 2 ,h] = Yl fi A hhh = (2p- l)!!A? iIaJs , 
{ii,...,i p }eii=i 

as claimed. □ 



Proof of Proposition 3.2. (a) The case #(/) = 4. This result can 
be established as a straightforward application of Lemmas 3.1-3.3. 

(b) The case = 6. Without loss of generality we can take / = {1, 2, 
. . . , 6} , and show that 

e ••• e nf h h h ) 

mn=— h m 6 3=-(3«=l 

E 5{ T MMM) = 0(li 1 ). 

7 Gr(6,3)\r P (6,3) 

It is sufficient to consider only the diagrams with no flat edges; it is readily 
seen that the latter can be partitioned into (a) the unconnected (unpaired) 
diagrams and (b) the connected diagrams. Now for (a) we note that, because 
there cannot be any flat edge, if the diagram is unconnected but not paired 
the set of rows must necessarily be partitioned into a group of two and a 
group of four; in other words, after some rearrangement of indices it must 
be possible to write any diagram 7 belonging to (a) as 7 = 71 U 72, where 
71 E Tp(2, 3) and 72 E Tc(4, 3). Thus, exactly as shown above, we obtain that 
the corresponding terms are bounded by C(2l\ + 1) . It suffices then to look 
at the connected diagrams. Consider first 7 E r c - i ( 2 )(6 ) 3); from Lemma 3.2 
and simple manipulations we obtain immediately 

(56) D[T CL{2) (6,3)]=0(l^D[T c (4,3)])=0(q 2 ). 
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Fig. 4. 7er a( 3)(6,3), 7 G F C£(4) (6,3). 



In case 7 G F c j^(6,3), it can be readily shown that we must have 7 € 

{Fcl(3) (6, 3) U r<7£(4)(6,3)}; in other words, these diagrams have at least a 
loop of order either 3 or 4 [F c j^r(6,3) is empty]. The latter claim can be 

established as follows. We have a graph with six vertices, each of which is of 
degree 3. Hence we can argue as in ([11], page 48) to show that the closure of 
the graph is complete (it is a clique). Then, from Theorem 4.5 of ([11], page 
48) it follows that the graph is Hamiltonian, that is, it contains a spanning 
cycle. With a permutation of the indices, we can then take the vertices to 
be at the corners of a regular hexagon, and we are left with an edge free 
for each of them. These edges will connect different vertices to form loops 
of order 3 or 4. The graphs corresponding to 7 G F cx( 3 )(6, 3) with at least a 
loop of order 3 are labeled type A; those with all loops of order at least 4 
[7 G r C £( 4 )(6,3)] are labeled type B (see Figure 4). 

Let 7 correspond to a type A graph; by Lemma 3.3 we easily have (see 
Figure 3) 

D[ r ,hMM = o(^l i f 2 ^)}^[r (7lM (4,3);Z 1) Z 2 ,Z3]) 

For type B, again up to a permutation of the indices we can take, with no 
loss of generality, 

mn = -m 2 i, mi2 = -77142, m 13 = -m 63 , 

m 22 = -m 32 , m 23 = -m 53 , m 3 i = -m 61 , 



m 33 = 



-m 43 , 77*41 = -77t 5 l, TO 52 = -m 62 . 



3G 



D. MARINUCCI 



The corresponding expected value is premultiplied by the factor 

(_l) m ll+ m 12+ m 13+m3i+m32+m33+m 5 i+m52+r?l53 _ J 



whence we have 

{ n Ae 

U=l,3,5 j=l my - 



(57) 



x 



( h h h 
\mu mi2 mi 3 

h h h \ f h h h 

niu -m 32 -m 53 J \m 3 i m 32 m 33 

h h h \ ( h h h 

m 5 i -mu -m 33 J \m 5 i m 52 m 53 

h h h 

m 3 i -m 52 -mi3 

^1 ^2 ^3 1 J h h h \ J ii ^2 ^3 



$:(-i)^(2x + i)P i 



3 ^3 x ) \h x h) \ x h h 

where we have used ([34], equation 10.2.3.17, page 339, and equation 10.2.4.20, 
page 340). In view of (55) we easily obtain 

(57) < -£= £ - = o(/r 1/2 / 2 " 1/2 / 3 1/2 io g / 3 ) = o^- 1 ). 

vnhh x= i x 

(c) The case #(/) = 8. Again, we can take I = {1, 2, . . . , 8} and show that 

h h 8 

e e n 



mu=— il m83=-^3 i=1 



^1 ^2 ^3 

mil m i2 m i3 



x £ <y( 7 ;/i,/2,/ 3 ) = oar 1 ). 

7 e{r(8,3)\r P (8,3)} 

As before, it is readily seen that the diagrams with no flat edges can be par- 
titioned into (a) the unconnected unpaired diagrams and (b) the connected 
diagrams. Now for (a) we note that 

7 G[{%(8,3)\r P (8,3)}nr F (8,3)] 

=> 7 g [{r c (/i, 3) r c (/ 2 , 3)} u {r c (/ 3 , 3) r c (h, 3)} 
u{r c (i 5 ,3)<g>r c (/ 6 ,3)®r c (j 7 ,3)}], 

where (Ii, 7 2 ), (I 3 , I4), are partitions of {1,2,..., 8} into disjoint 
sets such that 

#(/i) = 6, #(/ 2 ) = 2, 

#(/ 3 ) = 4, #(/ 4 )=4, 

#(/ 5 ) = 4, #(/ 6 ) = 2, #(/ 7 ) = 2; 
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by T\ <S> T 2 , we denote the set of all diagrams of the form 7 = 71 U 72, for 
71 G Ti and 72 G r 2 . For any two families of diagrams Ti(Ii, 3), r 2 (I 2 , 3) such 
that Ii H I2 = 0, it is readily checked that 

D[Ft r 2 ;h,l 2 ,l 3 ] = 0(D[T r ,h,l2,h} X D[T 2 ;h,l 2 ,l 3 }), 

D[r x u r 25 Zi, i 2 , z 3 ] = op[ri; Zi, z 2 , h] + £>[r 2 ; z x , z 2 , z 3 ] ); 

thus 

r>[%(8, 3) \r>(8, 3); l u l 2 ,h] 

= 0(D[r c (6,S)®r c (2,3);h,l 2 ,k}) 
+ 0(D[Tc(4,3)^r c (A,3)-h,l 2 ,l 3 }) 
+ 0(L>[r c (4, 3) r c (2, 3) (g> T c (2, 3); z lt z 2 , z 3 ]) 

= 0( J D[r c (6,3);Z 1 ,/ 2 ,y) 

+ 0(D 2 [r c (4, 3); Z 1; Z 2 , Z 3 ] ) + 0(D[T C (4, 3);h, l 2 ,h}) 
= 0(1^), 

as shown previously. It suffices then to look at the connected diagrams. 
Diagrams with a 2-loop can be handled by Lemma 3.2, and then by using 
results on lower-order diagrams. So we just have to consider 7 G ^^-^^(S, 3); 
these diagrams must have at least a loop of order 3 or 4 [in other words, 
r^-^y(8,3) is empty]. The latter claim can be established in the same man- 
ner as before; we repeat the argument for completeness. We have a graph 
with eight vertices, each of which is of degree 3. Hence we can argue as in 
([11], page 48) to show that the closure of the graph is complete. Then, from 
Theorem 4.5 of [11] it follows that the graph is Hamiltonian, that is, it con- 
tains a spanning cycle. With a permutation of the indices, we can then take 
the vertices to be at the corners of a regular octagon, and we are left with 
an edge free for each of them. These edges will connect different vertices 
to form loops of order 3, 4 or 5. Graphs with loops of order 3 can be dealt 
with as before through Lemma 3.3. If all loops are of even order, then the 
graph is bipartite ([11], Theorem 2.4, page 23); we label it a type C graph, 
for which we provide two isomorphic representations in Figure 5. 

If there is at least a loop of order 5, we label it a type D graph, for which 
two isomorphic representations are provided in Figure 6. 

To analyze the behavior of the components D[y, Zi, l 2 , Z 3 ] corresponding 
to type C graphs, we can take with no loss of generality 

m 11 = -m 21 , mi2 = -m 42 , m 13 = -m 83 , m 32 = -m 22 , 

m 3 i = -m 4 i, m 33 = -m 63 , m 51 =-m 61 , m 52 = -m 82 , 

m 53 = -m i3 , m 7 i = -m 81 , m 72 = -m 62 , m 73 = -m 23 , 
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Fig. 5. Isomorphic type C graphs. 




X 1 f Z2 



3 ^2 J I ^3 
^3 ^3 X \ ( l 3 k X 



h h h J I ^2 ^2 ^3 

where we have used ([34], equations 10.13.3.23 and 10.13.3.25, page 367); the 
sum runs over all positive integers x which satisfy the triangle inequalities 
h — h < x < li + 1\. By using (55), we obtain 

x 1 f h h x 1 f Z 3 Z 3 xl p 3 /i x 

^3 I2JW3 *1 k/U ^2 Z1JU2 ^2 *3 



(58) <C(2/ 3 + l) 2 max(i 2 i 1 

x [l 3 f 3 



In view of ([34], equations 10.13.1.1 and 10.13.1.3, page 361) the proof can 
be completed by an analogous argument for components corresponding to 
type D graphs. □ 
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Fig. 6. Isomorphic type D graphs. 
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